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Abstract. Gravitational radiation can be expressed in terms of an infinite series of 
radiative, symmetric trace-free (STF) multipole moments which can be connected to 
the behavior of the source. We consider a truncated model for gravitational radiation 
from binary systems in which each STF mass and current moment of order I is given 
by the lowest-order, Newtonian-like /-pole moment of the orbiting masses; we neglect 
post-Newtonian corrections to each STF moment. Specializing to orbits which are 
circular (apart from the radiation-induced inspiral), we find an explicit infinite series 
for the energy flux in powers of v/c, where v is the orbital velocity. We show that 
the series converges for all values v/c < 2/e when one mass is much smaller than 
the other, and v/c < 4/e for equal masses, where e is the base of natural logarithms. 
These values include all physically relevant values for compact binary inspiral. This 
convergence cannot indicate whether or not the full series (obtained from the exact 
moments) will converge. But if the full series does not converge, our analysis shows 
that this failure to converge does not originate from summing over the Newtonian 
part of the multipole moments. 
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1. Introduction 

The possibility of detection of gravitational waves from inspiralling compact binaries 
using laser interferometric gravitational-wave observatories such as the U.S. LIGO 
and the French-Italian VIRGO projects has brought into sharp focus the accuracy 
of calculations of gravitational waves using approximation methods. The ability 
to measure the source parameters using matched filtering of theoretical templates 
against the tens of thousands of cycles observed in, say, a double neutron-star inspiral, 
assumes that the templates are sufficiently accurate, especially in the evolution of the 
phase of the waves, that the errors are smaller than the errors arising from noise 
in the detectors |). This may require knowledge of the damping of the orbit 
via gravitational-radiation reaction (which determines the non-linear evolution of the 
phase) that incorporates corrections to the lowest-order quadrupole approximation as 
high as order (v/c) 6 || |[. Corrections at 0[(v/c) 4 ] and 0[(v/c) 5 } have recently been 
calculated || ||, 0, ||. 

On the other hand, there is evidence that such post-Newtonian expansions of 
gravitational radiation (weak-field, slow-motion expansions in powers of e 1 / 2 ~ (v/c) ~ 
(Gm/rc 2 ) 1 / 2 ) do not converge rapidly, if at all. For the case of a test body in circular 
orbit around a black hole, perturbation calculations carried to very high order in v/c 
show slow convergence — the coefficients of successive powers of v/c in the expansions 
grow alarmingly quickly |l(J |LT|] . This could call into question the accuracy of any 
post-Newtonian approximation truncated at a finite order, and by implication, the 
use of such approximations in templates used in data analysis for LIGO and VIRGO. 

To date, no method has been identified to study the convergence properties of 
the post-Newtonian expansion rigorously and in generality, because of the difficulty 
in generating higher-order corrections explicitly. In the case of post-Newtonian 
calculations for systems of arbitrary masses || [?J , the complexity of the computations 
grows rapidly with each succeeding order. Although the "third" post-Newtonian 
order (3PN), corresponding to corrections at 0[(v/c) 6 ] may be achievable, progress 
beyond that is unlikely. Besides, 3PN order may be adequate from the data-analysis 
point of view, given the expected noise characteristics of the advanced LIGO/VIRGO 
detectors, provided one had some understanding of the convergence properties, or 
could bound the errors neglected. 

There is one situation in which the convergence properties can be studied at 
all orders in v/c. We call this the "bare-multipole truncation", and although it is 
unphysical, it may provide clues concerning the convergence of the true (and physical) 
post-Newtonian expansion. This truncation is best discussed with reference to figure 
1. 

Gravitational radiation can be expressed in terms of an infinite set of radiative 
mass and current "symmetric, trace-free (STF)" multipole moments, which can be 
related to multipole moments of the source |Q . For example, the rate of energy loss 
E and the gravitational waveform far from the source /i^ T can be written 

^=^E( a; « + Vl ai " a '« + Vl ai " a ' +ft C+i)j ai - Q! ('+i)j Ql - ai ) , (la) 

C 1=2 

**r - £?t (J - + ^ e „ ( , ■) (1« 
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with 

(Z + !)(/ + 2) 4Z (Z + 2) 

1 I(J-l)Zl(2J + l)!!' n (i-l)(Z + 1)1(21 + 1)!! ' M 

where and J ai --- a ' are respectively the STF mass and current multipolc 

moments of order I, 7V ai ' °' denotes an /-dimensional product of unit vectors pointing 
from the source's center of mass to the observer at a distance R, e pq j is the completely 
antisymmetric Levi-Civita symbol, and the subscript "TT" denotes the transverse- 
traceless part. The indices on the left of the moments denote a number of derivatives 
with respect to retarded time. Summation over repeated indices is implied. 

Each STF multipole moment is then expressed as a post-Newtonian (PN) 
expansion, effectively in powers of e 1 / 2 , starting with a lowest-order multipole moment 
given in general by 




where tha, x\, and v % A are the suitably defined mass, position, and velocity of each 
body, and the superscript notation < I > is short-hand for the I indices. 

Figure 1 represents these STF moments in an array with multipole order I 
increasing horizontally beginning with I = 2 (for each /, I <1> and J <(_1> are grouped 
together), and with PN order for each multipole moment increasing downward. (Note 
that I <0> , 1 <1> , and J <x> are related to the non-radiative total mass, center of 
mass, and total angular momentum of the system.) It is straightforward to show that 
the first-order PN correction to each multipole moment is O(e), followed by a "tail" 
correction at 0(e 3 / 2 ), then 0(e 2 ), and so on. We ignore finite-size effects, such as 
the effects of spin, in this discussion. There is evidence from black-hole perturbation 
studies that the PN expansions become non-analytic at high enough order, with the 
appearance of In e terms (|o| [ll] . 

Because each multipole moment is differentiated with respect to retarded time a 
number of times corresponding to its multipole order, and because d/dt ~ e 1 / 2 /^, a 
specific moment contributes to hl£ T and to E at a PN order related to its multipole 
order. For example, a determination of E to 2PN order requires knowledge only of 
the first three columns of figure 1, (/ = 2, 3, 4), together with corrections through 
2PN order of the quadrupole STF mass moment (1 <2> ), the 1PN corrections of the 
octopole mass and quadrupole current moments (I <3> , J <2> ), and the lowest-order 
hexadecapole mass moment (I <4> ) and octopole current moment (J <3> ). These 
correspond to the shaded region in figure 1. For the waveform at 2PN order, on 
the other hand, the moments and their corrections indicated by the dark border are 
required. The difference derives from the fact that the waveform is linear in the 
moments, while the energy flux is quadratic. 

The "bare-multipole truncation" consists in keeping only the first row of figure 1, 
but keeping all orders in I. This is not a consistent PN approximation to cither the 
waveform or the energy loss. Moreover, this truncation does not correspond to any 
sort of physical model for the system, such as one in which non-gravitational forces are 
responsible for keeping the two bodies in orbital motion. This can be seen as follows: 
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Consider the O(e) corrections to a particular multipole moment. These contain terms 
of the form Gm/rc 2 , which can be thought of as the gravitational contribution to the 
effective stress-energy tensor, and terms of the form (v/c) 2 , which can be thought of as 
special relativistic corrections in the matter contribution to the effective stress-energy 
tensor. Now, it is clear that the first set of correction terms would be discarded in 
a non-gravitational model for the source motion, although it would presumably be 
replaced by something else (related the stresses responsible for maintaining the orbit). 
However, it is also clear that the second set of terms cannot be discarded in any sort of 
physical model for the source. The fact that both sets arc discarded in our truncation 
clearly makes such a procedure unphysical. 

Although it is unphysical, the bare-multipole truncation may illustrate some of 
the convergence properties of the full formulae. For two-body systems, we show that 
the bare moments have the general form 



J<f =^vf l+1 (r)) LoS'" 1 _ , (4) 

J Tr 



where /i = mim2/(mi + 7712) is the reduced mass, 77 = ///(mi + ma), Sq is a 
symmetrized product of I unit radial vectors directed from body 2 to body 1, L is 
a unit vector in the direction of the orbital angular momentum, r and v are the 
magnitudes of the relative separation and orbital velocity, respectively, and o denotes 
a symmetrized product. The function fi(j]) is given by 

m = [ P i - i +(-i) 1 ] /(1 + P ) i -\ (5) 



p = m 2 /m 1 = — 1-277-^/1-477 , (6) 
Z77 L J 

where we choose the convention 777,2 < mi, so that < p < 1. 

Specializing to quasi-circular orbits with angular velocity f2 = v/r, we find, after 
considerable manipulations of STF tensors (Section 3), the closed- form result 

E = E Q \l + T i B l f 2 ( V )(l) 2l ' i ) , (7) 



(=3 



where Eq represents the quadrupole approximation, 

£ Q = f^W. (8) 

The coefficients Bi are given explicitly by (p8|) . It is interesting to note here that (7) 
does not depend on any assumption of an equation of motion for the bodies; only the 
existence of a quasi-circular orbit with angular velocity Q is assumed. 

A completely equivalent series, applicable in the test mass limit (77 — * 0, |/z(t7)| — ► 
1) may be derived using black- hole perturbation theory (see Section 4 for details). 
The coefficients Bi resulting from that method have a different, albeit numerically 
equivalent representation, given by (|44|, [|5|). Most importantly, they can easily be 
shown to have the form, for large I, 

5 f l\ 1/2 /ex 2 ' 



s -s ; (I) MM- < 9 > 
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where e is the base of natural logarithms. 

By applying the standard Cauchy ratio test, we find that the radius of convergence 
of the sequence is given by 

/ v \ 4 

Vc/convorgo c(l + ^1- 4r/) 

The series thus converges for values of v/c less than a critical value ranging from 0.74 
{rj = 0) to 1.47 (i] = 1/4). These values encompass all physically relevant values in 
binary inspiral, until the final coalescence phase. 

The infinite series (7) can be used to study the error made in truncating the 
method at a given PN order. (By this we mean a truncation of the sum appearing 
in Eq. (7) at the value of I corresponding to the specified PN order. For example, a 
truncation at 2PN order involves keeping only the I = 3 and I = 4 terms.) Figure 2 
shows, as a function of v, the fractional difference between a given PN truncation and 
the full series. Notice that, at v ~ 0.4, corresponding to the innermost stable orbit 
for a test-body orbiting a black hole, the errors at 2PN and 3PN order are 2.1 and 
0.6 percent, respectively. For equal masses, odd-numbered mass moments and even- 
number current moments vanish, thus there are no terms in E at odd-PN order in the 
bare-multipolc truncation. This is partly responsible for the improved convergence 
in the equal-mass case: at v ~ 0.4, the fractional difference between the 2PN/3PN 
approximation and the full series is 2 x 10~ 5 . 

The remainder of this paper provides details of the calculations. In Section 2, 
we describe the bare-multipole truncation and set up our definitions and conventions. 
Section 3 derives the energy flux, and discusses the convergence of the series. In Section 
4, we describe the analogous approach using black hole perturbation theory, and 
find an analytic expression for the radius of convergence. Section 5 gives concluding 
remarks. A number of technical details are relegated to Appendices. 



2. Bare-multipole truncation for gravitational radiation 

We describe the motion of a system of two bodies whose size is negligible compared to 
their relative separation. We work in a coordinate system whose origin is the center 
of mass, and define the relative position vector x = xi — X2, where xi and X2 are 
the positions of each body. If the masses of the bodies are m\ and 1112, we define 
m = mi + JTt2, \i = mi m2/m, rj = n/m, and dm = mi — ni2. Also, r = |x| and v = x. 
Henceforth, we use units in which G = c = 1. 

The radiative multipole moments in (la) and (16) are expanded in powers of 
e 1 / 2 ~ v <~ (m/r) 1 / 2 ~ r, in the generic form 

I< 1 > ~I<j>(l + e + e 3 / 2 + e 2 + ...), (11) 

with a similar expansion for 3 <l> , where the e 3 / 2 contribution signals the first 
appearance of the effects of radiative "tails" , caused by backscatter of the waves off the 
background spacetime curvature. In calculating E, one takes each multipole moment 
of order I, differentiates it I + 1 times with respect to retarded time, and squares the 
result (contracting on all indices). Because the /-pole moment is differentiated once 
more than the (I — l)-pole moment, its leading contribution to E is of order v 2 higher 
than that of the (Z — l)-pole moment. Additionally, the current multipole moments 
are proportional to the angular momentum of the system, through a factor £j a & x a v b , 
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which adds an extra velocity v b to their expressions; consequently the contribution of 
a current multipole moment is 0{v 2 ) higher than that of the mass multipole moment 
of the same rank. As a result, (la) is a post-Newtonian expansion of the luminosity 
in powers of e. Similar considerations apply to the expressions for the gravitational 
waveform. 

Explicit expressions for the radiative mass and current multipole moments for 
general two-body systems are now known sufficiently accurately to calculate the energy 
flux to 5/2PN order beyond the quadrupole approximation ||, and the waveform to 
2PN order || . By way of illustration, we quote here the moments sufficient to calculate 
the 2PN energy flux §, §: 



/'< 



+x l x 3 



l + i(l-3^-i(5-8,)(^) 



- — (1 — 37/) r r x l v J + — (1 — 3rj) r 2 v* v 3 



504 
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(2021 - 5947 7/ - 4883 rf 
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(742 - 335 
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+ r Tail ) 
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x l x 3 x k 
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£iab I 1 + \ (1 - 577) v 2 + 2 (1 + 7/) (™) 



STF 



rra- 7 ) x a u b l 



STF 



rijfe _ 



= (1-37/) (e mba; a w Vx fe ) 



STF ■ 



(12) 

Explicit expressions for the "Tail" term may be found in (^,0], for example. 

We now specialize to quasi-circular orbits, i.e. orbits which are circular, apart 
from the slow inspiral caused by radiation damping (at sufficiently high PN order the 
non-circularity of the orbits must be taken into account). We thus approximate: r = 0, 
x = 7'fi, v = v\, where n and A are unit vectors in the radial and tangential directions, 
respectively, and nxA = L. We define the following symmetric tensors (see Appendix 



A): 



. X ik h ik+1 h lk + 2 . . 



(13) 



where the sum is over all possible distinct permutations Af 



of the set of 



indices {ij}. The STF tensors I <(> and J <z> can then be expressed in terms 
of these symmetric tensors, and the post-Newtonian terms can be converted, via 
the appropriate equations of motion, into terms involving v alone. The result is, 
schematically, 



j<2> 
j<3> 
j<4> 
j<2> 

j<3> 



M r 2 [(l 

3 



OtiV 



a 2 v A )S 2 +v 2 {l + a 3 v 2 )S° 2 ] 



STF 



, T<2> 
+ r Tail ' 



= fir 



[(1 + ^)8 



rn 

M 7 " 4 [7iS 4 ] STF 
2 / (5m 



/3 2 f 2 S^] STF 



LoS* 
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LoS; 
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(14) 



where the coefficients a*, /3j, etc., result from combining expressions in (12), and 
depend only on 77; o denotes a symmetrized product, defined by 



LoS 



T-l _ 



EE L^ff 



1 



\ " 



T ^ 



(15) 



where ILp are all the T distinct permutations of the indices {ij}. We now restrict 
these expressions to the leading order in v, and denote them the "bare" multipole 
moments 1 5,^ and J 5,^ , as follows 



T <2> 

*(0) 

T <3> 

*(0) 

T <4> 

*(0) 

T<2> 

J (0) 



I<3> 
rf (0) 
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STF 
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STF 
STF ' 



fir 2 v 



LoSj 



STF 



fir 2 vf 3 (rj) 



to Si 



STF 



/ir 3 i>(l — 3rj) 



LoS 2 , 



— ,,^,3 



STF 



fir u/ 4 (r?) 



LoS 2 , 



STF 



(16) 



where the function /;(ri) is given by (6) and is plotted in figure 3 for different values of 
j] within the possible range (0,0.25]. No te that \fi(rj) \ < 1, for all 77 7^ 0; |/;(0)| = 1, 
for all I and /;(0.25) = for odd I (see Appendix B| ). Equations (16) suggest, and 

^> and T<*> 
(o) ana J (0) 



Appendix B confirms, that the general forms for tj^F* and ifS? at leading order are 



given by (4) . Notice that the original expressions for 1 5F and J ^F were expressed 
in harmonic coordinates. However, these expressions are also valid in Schwarzschild 
coordinates, as the transformation between coordinate systems, both for r and v, 
introduces terms that contribute only higher-order corrections. 



Finally, we explicitly extract the traces in the previous expressions, and get f(o) > 
and J(^ ) ) > in terms of the symmetric tensors S Z fc _fe only, 



<i> 

(0) 



F'/,WE ^P (D W <^ (17«) 



Jg> =Mr',/mW t ( ' X) ^ Q ( ^ 2J) (D ( ^(LoS»H)) . (176) 



Here [a/6J denotes the integer part of a/6, and is a symmetrized product of j 
Kronecker deltas, given by 

jyU) — gih'hfrh'U ... J*2i-l,*2i) _ 1 ( J»l»2j»3»4 ... ^*2i-l,*2j ^8) 

■A' (ii 2) n 



where 



j\\2j\ 2 J \2 

is the number of distinct products of deltas. The factor (I — 2j)/l in (176) comes 
from taking into account the fact that the unit vector L is perpendicular to n when 
applying traces to the tensor L o S _2j ' _1 . 

Notice that each term in (17a) and (176) has rank I — 2j, varying for each j in 
the sum. This fact complicates considerably the process of multiplying two of these 
objects, as prescribed in (la). In the next section we will describe how to simplify 
this calculation. 



3. Energy flux in the bare- mult ipole truncation 

We now express the energy flux approximately in terms of the bare multipole moments 



oc 



<i> (*+i)t<;> 
(o) • •>) 



(20) 



where the notation (•) stands for the inner product that saturates all the indices 
{ii} = {ii, i 2 , ■ ■ ■ , ii} of the two tensors. As we are dealing henceforth with bare 
multipole moments, we drop the subscript (0) notation. 

In order to simplify the product of derivatives in each of the sums of (20), we first 
write: 

(i+i) t <;> . (i+i)p<;> _ ^_ [(Of <«> . (j+i)j<i>1 _ (0f<«> . (i+2)j<;> 
dtl \ 



d_ 

It 

d_ 

It 

d_ 

It 



d_ 
~dt 



(i-i)f </> . ((+2) j<;> 

] + ( i - 1 )f<'> . (*+3)j<«> 

] + (-1)' +1 I < ' > • (2'+ 2 )f <'> 



(i-l)f <l> . (*+3)j <l> 



(21) 
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The first term is the total derivative of several terms, each of which is the product 
of an even and an odd number of derivatives of 1 <1> , respectively. Since each time 
derivative converts n into A, each term is thus a product of a tensor with an even 
number of A unit vectors by a tensor with an odd number of A's; consequently, every 
term in the bracket is null, because n- A = 0. The same can be done with the magnetic 
terms, so that we may write 



(i + l)j<i> . (Z+l)j<;> _ (-l) / + 1 J<'> . (2/+2)j<i> 
(i+l)j<i> . (l+l)j<l> _ 



i+i i</> . (2i+2) t<;> 



(22) 



Now, 1 <l> and 3 <l> are trace-free tensors, as are their derivatives. But the 
terms in \ <l> that contain Kronecker deltas, such as 5 %1%2 (Rest) 13 '" H , are just trace 
operators on the subspace {ii,Z2}; so their product with the time-differentiated STF 
tensor ( 2i + 2 )f< ! > w in vanish. Thus, the only non-vanishing term is the symmetric 
product of / unit vectors [h 11 h %2 ■ ■ ■ multiplied by ( 2i + 2 )f <'>. A similar argument 
applies to the current moments. The result is 

(i+i) t <;> . (i+i)j<;> = (-ly+i^ [^1^2 . ..^tjj . (2i+2)j<i> 



(-l) l+1 fir l S l ^ 2l+2 h <l> 



(i+i)j<i> . (i+i)j<;> 



(-l) l+1 fj,r l v 



. (2Z+2) j<Z> 



= (-1 



) l+1 (j,r l v 



L o S^ 1 



, (2Z+2) j<i> 



(23) 



We now need an expression for the (21 + 2)-time-derivative of l <l> and J <z> . 
These tensors (17a, 176) are sums of products of constants with respect to time (like 
D( J ) and L) with a symmetric tensor Sq 2 "' , which has only radial components. We 
will concentrate first on calculating the 2p = 21 + 2 derivatives of a symmetric tensor 

where T = I - 2j = I, I - 2 ,1 - 4, • • •. 

Th e derivative of any symmetric product of unit vectors, radial and/or tangential 
is (see Appendix A.l) 



= n [(i - k)s 



k+l 



(24) 



Thus, even if we have only radial components in Sj, tangential components will appear 
after differentiation. Notice that ( g4j ) is the same rule as the derivative of a product 
of k cosine and I — k sine functions 



dx 



(sin' k ilx cos k Ox) = 



(/ - k) (sin'- fc - 1 fix cos k+1 fix) 



k (si 



fix cos*" 1 fix) 



(25) 



We may thus simplify the calculation considerably by associating each symmetric 
tensor with a product of sines and cosines in the form 



sin' k fix cos k fix. 



(26) 



Any higher-order derivative of S k , expanded in terms of others S fc , , will have 
the same coefficients as the corresponding derivative of the analogous product of sines 
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and cosines, expanded in products of sines and cosines, provided that we can express 



such a derivative in terms of the equivalent products of sines and cosines (in Appendix 
~] we show that there is a subtlety in this procedure) . 
We are interested in the (2p)-derivative of Sj, which we associate with sin T fix, 
with T = I, I — 2, I — 4, • • •. To get a closed expression, we first expand sin T fix as a 
sum of sines of multiple angles (l3) : 

{[___J "I 
E L¥J ~ k (I) cos ~ 2k > + \ ( j ) \ > if T = 2n, 

1 I T-l I /T\ 

sin T fix = - frT ^ (-l)L— J+M j sinfi(T- 2fc)x, ifT = 2n-l. (27) 

We then take the 2p derivative of (27): 

^ (sin T fix) = ( ~ 2 ^ 2P 2 L * J ~* (T) ( T - 2fc ) 2P cos " 2k ) x ' if T = 2n ' 

fc=0 ^ ' 

^(sm T Vlx)= £ (-1)^+" (T) (T-2fc) 2p s infi(T-2fc)a:, 

if T = 2n - (28) 

Now we can express cosi7(T — 2k)x and smQ(T — 2k)x back in terms of powers of 
sines and cosines, to recover terms analogous to the symmetric products Sjr fe : 

cosP'" 



ifl(T-2k)x = E C 2 q k ) sin2?fix cos (T - 2fe - 29) fix, \iT = 2n, 



q =o ^ 



sinfl{T-2k)x = sin fix E (-!)« ( g + j ) sin 29 fix cos (T - 2fc_29_1) 



9 =o 



if T = 2n {29) 

Combining (28) and (29), one can get the 2p-derivative of a product of T sines. The 
2p-derivative of Sq will then have a similar expansion, with identical coefficients: 

(2p) (SI) - E J (-D Li J "* (I) ( T - 2 *) 2P 

x E 2 , ) S T- 2fc - 2? > ifT = 2n, (30a) 

[___.J 

(2p) (S^) = ^g^ E (-l)L^WT)(T-2*)* 
fc=0 ^ ' 

x 



9=0 



E (-^'Ug + iJS^Wi' if^ = 2n-l. (306) 
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However, this expansion of (2p) (So) is not yet what we need. Notice that 



preserves the rank of the tensors involved: rank(S' fc ) 



the derivative rule (\24 
= ran 

2k, different for each k. This is an artifact created by expanding the 



rank(S / fe+ ' t 1 1 ) = rank(S| c _ A j +1 ). This is no longer true for (30a, 306), where each term 



is of rank T 

sin T £lx in terms of sines and cosines of multiple angles, which we resolve in lAppendix 



The product of these expressions with S or with Lo Sq 1 in (23), selects only the 
terms of the form Sq in equations (30 a, 306). The expressions simplify considerably 



(see Appendix A. 3), giving the same result for even or odd T. Using (A23), we get 
the following product, valid for any T, T' (even or odd) and for any p: 



where 



E 

k=0 



(T - 2k) 2p . 



(31) 



(32) 



As the angular momentum L or the Kronecker delta operators are constants with 
respect to time derivatives, we also find that 



aT' 
°0 



A(T,p) (Sjf 



D«) o S^ 



LoS 



= A(T,p) 



LoS 



T' 



D 0') 



LoS 



(33) 



With these simplifications we can now express the luminosity in terms of products 
of symmetric tensors containing only radial components. This makes it possible to 
calculate these products in closed form, for any multip ole order, I. Using the previous 
expressions and the following rules (see Appendix A. 2): 



DW o S " 2j 



1, 



LoS 



l-i 



D 0) 



LoS 



the products (23), with T = I — 2j, 2p — 21 + 2 and v = fir are 

(;+Dj<;> . (*+i)j<*> = 2- ( - l - 1 ^ 2 v 2l n 2 tf(r 1 )C l , 
(i+1)j<;> _ ( i+ l) J<; > = 2 -(i- 2 ) r 2 A1 2 u 2z+2 n 2 f? (rfi-Dt 



where C\ and T>\ are given by 



W (-i) J ^ C-) Q ^ 

l^i (2l\ 2—1 

j=0 \2j) k=0 



l-2j 



V, 



E tl^Ml (! _ 2 ,, 



E 

k=0 



- 2j - 2k) 2l+2 
l-\-2j 



(34) 



(35) 



(36a) 



3=0 

Finally, the luminosity in terms of Ci, T>i, and fi(rj) is 



(I - 1 - 2j - 2fc) 2 '+f366) 



E : 



1 



1=3 



Bif?(v)v 



21-4 



(37) 
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where Eq 

by 



-/it 2 t> 4 f2 2 is the luminosity due to the quadrupole term and Bi is given 



B, 



5 (I + 1)(Z + 2) 
16 Z(Z- 1) (2Z + 1)! 



16(2Z + 1)Z 

(Z-2) (Z + 2) 



(38) 



To study the convergence of this series, we apply the standard Cauchy ratio test, 
requiring that 



urn — v < 1 

i-oo Biff{rj) 



(39) 



Evaluating this numerically up to I = 250, we find the constraint \fi+i(j])/ fi(r])\ v — ► 
(1 + V 1 ~ 4r?)u/2 < 0.74. The case r\ — 0.25 must be treated separately, since 
J; (0.25) = for odd Z. For this case, the Cauchy ratio test takes the form 
\fi+2/ fi] 1 ^ 2 v — ► v/2 < 0.74. which is the continuous limit of the previous statement. 
Note that, for a given v, the convergence is worst for rj = (test mass limit) and best 
for 77 = 1/4 (equal masses). 



4. Black-hole perturbation theory and the bare-multipole truncation 

In this section we consider the calculation of the bare multipole moments in the 
restricted context of a binary system with small mass ratio. We shall therefore demand 
77 <gC 1 throughout this section. The method of calculation used here is completely 
different from the one used in the previous two sections. 

For the problem considered in this paper, the internal structure of the orbiting 
masses is of no consequence. For simplicity, in this section we take the larger mass m\ 
to be a nonrotating black hole. The smaller mass m,2 then creates a small perturbation 
in the gravitational field of the black hole. This perturbation propagates away from 
the source as a gravitational wave. By virtue of the small mass approximation, the 
gravitational perturbations are accurately described by solving a linear wave equation 
on the Schwarzschild spacetime. This equation is called the Teukolsky equation fl4| , 
and it can be solved exactly, for example, using numerical methods |To| . 

We assume that the mass mi moves on a circular orbit with radius r, such that 
v = (m/r) 1 / 2 is much smaller than unity. (Here, m = mi + mi — m\ is the total 
mass.) In this limit the Teukolsky equation can be solved analytically, and the bare 
multipole moments of the radiative field can be evaluated. Such a calculation was first 
carried out by Poisson |ll| , and then completed by Poisson and Sasaki |l(| , hereafter 
referred to as PS. 

The Teukolsky equation is analyzed by first separating the variables. This means 
that the radiative field is expressed not in terms of symmetric trace-free tensors, but 
in terms of (spin- weighted) spherical harmonics ]l7| . The two representations are 
entirely equivalent |M , and in both cases the gravitational- wave luminosity takes the 
form of (7). Here, of course, rj — > so that \fi(r])\ = 1. 

The bare multipole moments were calculated explicitly by PS, who express the 
gravitational-wave luminosity as 

£ = ^q i+EE^l (40) 

\ 1=3 m=l ' 
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apart from a slight change of notation; cf. their equation (5.14). The contributions 
rji m take different forms according to whether I + m is odd or even: 

1 for / + m even 



nil. 



Plm V 
qirn V 



21-2 



for I + m odd 



(41) 



The coefficients pi m and q\ m are explicitly given in equations (5.15) and (5.16) of PS. 
Simple manipulations bring (J4C|) into the form of (7), which we write as 



E = Eq ( 1 



5> 

1=3 



.21-4 



(42) 



We find 



JY 



Bl = Pll + y~] (Plm + Ql-l,m) ■ ( 43 ) 
n=0 

Here, m is to be considered to be a function of n: for odd I, m — 2n + 1 and 
N = (I — 3)/2; for even I, m — 2n + 2 and N = (I — 4)/2. This ensures that the sum 
is properly restricted to values of m such that I + m is even. The largest value of m 
contributing to the sum is therefore 1 — 2. 

Some algebra, using the explicit expressions for pi m and qi m provided by PS, 
converts ( fi"3| ) into 

5(Z + !)(/ + 2) 



Bi 



16(2/ + l)!(/-l)i 



21+2 



where 



n 2 



N 



m M+2 {l-m)\ (I + m)\ 
{2l)\l 2l + 2 4^ [(l-m)/2]\ 2 [{l + m)/2}\ 



1 



4(2/ - l)(2l + - m)(l + m) 
(l-2)(l + 2)m 2 



(44) 



.(45) 



This expression for Bi is entirely equivalent to that given in (|38|). We have indeed 
verified that these Bi's are numerically equal to the ones derived in section 3. However, 
we have not been able to establish their equality algebraically. 

We now wish to examine the convergence of the sequence B\v 21 . This analysis 
will greatly benefit from the simplicity of our current expression for Bi . We point out 
first that the sequence 6/ converges. This was established by numerical experiment, 
which also reveals that b^ ~ 0.01. This implies that the behavior of Bi as I — > oo 
is determined entirely by the factor l 2l+2 /(2l + 1)! in ([44). Writing (21 + 1)! as 
(2/) 2 r(2Z)[l + 0(1/0] and using the Stirling approximation 18 

(2^) 1 / 2 e" 



T(z) = (2TT) 1 ^e- z z z - 1 / 2 [l- 
we quickly arrive at the asymptotic form 



0(V*)] , 



Bi 



5(l + froo ) (I 

64 



1/2 



i + o(r 



(46) 



(47) 



We therefore find that the sequence Biv 21 behaves asymptotically as (ev/2) 21 . A 
direct application of the Cauchy ratio test then reveals that the sequence converges 
provided 

v < 2/e ~ 0.7358. (48) 

The generalization of this result, appropriate for the case of nonvanishing mass ratios, 
was given in section 1. 



14 



5. Concluding remarks 

We have shown that a truncated model for gravitational radiation from binary systems 
in circular orbits converges for values of the orbital velocity that encompass all inspirals 
of physical interest. However, our model is admittedly non-physical, and may only 
be revealing that, whatever poor convergence properties have been seen to date in 
the PN expansion, they do not arise from summing over the Newtonian part of the 
multipole moments, but arise instead from the PN corrections to the moments. 

Figure 4 illustrates the limitations of our truncation and reinforces the notion that 
it is unphysical. Shown is a comparison, in the test body limit, between our bare- 
multipole series and the physical results of black-hole perturbation theory, including 
both the numerical "exact" results Q] and the true PN series, accurate to 3.5PN order 
Pf . It is apparent that the truncated model compares rather poorly; this is mostly 
due to the fact that while our series is a sum of positive terms, the true PN series 
is alternating. Nevertheless, our analysis can be improved somewhat. The lowest- 
order tail corrections to each STF moment can be written down explicitly, and could 
thus be added in a straightforward way. It is possible that the O(e) corrections to 
each moment could also be calculated without too much difficulty. Whether such a 
"dressed-multipole" truncation shows better agreement with the exact results — and 
still converges — will be a subject for future work. 

Appendix A. Symmetric products of unit vectors S l ^ k 

The multipole moments in (12) can be expressed in terms of products of components 
of unit vectors, 

(Pi - *)* 1 ""*' = A* 1 A 12 • • • X ik n ik+1 n ik+2 ■ ■ ■ h il . (Al) 

Here k components are in the tangential direction to the circular orbit, A*, and (I — k) 
are in the radial direction, n l , with nn = n % h l = 1, A- A = A 1 A 1 = 1, n- A = n l X l = 0, 
and n x A = L (or L l = eijkW\ k ). 

A symmetric tensor can be constructed from (Pj i T' c ) 1 



(Sfc - *)* 1 '"* 1 = " Xin " y2 ■■■X lk n lk+1 n lk+2 ■■■h ll) 
1 

= Jf 



\x il X 12 ■ ■■X lk h lk + 1 h lk + 2 ■■■h 11 



= ^E( p *"*) <1 '" i ' ' ( A2 ) 

where the sum is over all the M = (,) permutations of the indices {«/}. 

For simplicity, from now on we will drop the set of indices {i± ■••«;} in the notation 
for P l k ~ k and S l k ~ k . 

Appendix A.l. Differentiation rule for Sj,r fe 

In a system of two bodies following a circular orbit, the time derivative of the 
components of the normal and tangential unit vectors are, respectively, 
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where f2 is the angular velocity, such that v = fir. 

We need to calculate the time derivative of a symmetrized product of I components 
of unit vectors, where k of them are tangential (As) and (I — k) are normal (h l, s), as 
in (A2). The time derivative of one of the Pjr fc 's that belongs to Sj.r fe : 



dt ^ k 



l-k\ 



fi 



n %1 X 12 ■ ■ ■ X th n lk+1 h lk + 2 ■■■n 1 ' 
+ \ n rY 2 ■ ■ ■ X ik n Lk+1 h lk + 2 ■■■n 1 ' 



+ \ n '\' 2 ■■■h lk h lk +^h lk + 2 ■ ■ -h il \ 
+ (\ H X 12 ■■■X k X^ 1 h lk + 2 ■ ■ ■ h H 
+ \ ll X 2 ■■■X k h lk+1 X lk + 2 ■■■n il 



+ X 11 X 12 ---X k h ik + 1 n ik + a •••A* ! ) 



(A4) 



is the sum of two parts: the first parenthesis with k terms of the type PjfLj , and 
the second parenthesis with (/ — k) terms of the type P 

Taking the derivative of the whole Sjr fc creates (?) sets of terms like (A4). In the 



l-k-l 
k+1 ■ 



expression for d (Si ) /dt, all permutations of indices will be present because they 
were present in the original S l k 7 k - We just need to collect all the terms of the type 



k 

k _ 1 , that belong to S^l^ +i , and all the terms of the type P^p 1 , that belong to 



There are ( fe l 1 ) P 1 fe _* +i -terms and ( fc+1 ) P^+i 1 -terms. So the expression 



for the time derivative of Si k is 



(fc+i) 



n [(i - k) stff 1 



• pi-fe-i 
-ks[-_\ +1 ]. 



(fe-l) 11, 



/ j r fe-l 
{ij} 



(A5) 



For this procedure to be valid, at least one of each of the new symmetrized products 
\ +1 has t 

)(l-k) 



S L k+i 1 an d ^fc- fe i +1 nas to be obtained. In other words, 



(fc+i) 



> 1, 



and 



k 



> i . 



(A6) 



which are trivially true, with the equality corresponding to k = I (for S;), and k = 
(for Sq), the two possible limit cases. Then, (A5) is the rule of differentiation for the 
symmetric tensors S l ^ k . 

Appendix A. 2. Products of Sq 's 

In this section we will limit the discussion to symmetric products that contain only 
radial components. In this case, (with the same conventions as in (|l5|), (18), and (19)), 
we have 



ST = n (tl n' 2 ■ ■ -h lT) 



(A7) 
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LoS 



T-l 



D (i) o Sq~ 2: > = 6 {lll2 5 l3U ■ ..5 i2 *- 1 > i2 i'n i2 *+ 1 ■ ■ -h lT) 



AA(T;i,2; 



(A8) 



(A9) 



L O S n 



^(ii^jiau . . . gi2j-i,i2j fisj+iffiaj+s . . . fair) 
1 



(T-2j)JV(T-l;j,2) 



V* 5 hi2 ■ ■ • <J*«-i '*y£«i+i n i2 ^+ 2 • • • n* T , (A10) 



where II in the latter two equations simply denotes the distinct permutations, and 

(All) 



T-2j 



Using the orthonormality of n and L, and the definition ■ nn = 5 % i n % n> = 1, we 



L O S n 



Lo S 



T-l 




calculate the products: [DW o S T 2j ] • Sg and 

which occur in (23). Notice that we use as "target" (or second factor in the product), 
the simpler of the two factors. This simplifies the counting of resulting terms. For 
both cases, the expression is, schematically, 



— ^Terrm • -L^Term 2 



(A12) 



where jVj is the number of distinct terms in each sum. Next, we calculate the effect of 
multiplying one term only from the symmetric first factor by the whole of the target, 
say 



[Termi] 



There are A/i identical products like this in the final result, so 



-.^Terni! • -L^Term 2 



[Termi] 



[Afi ^ J [AT 2 

For the first product, flA13j ) gives 

trivially, as the target has only one term. The result is 



7-^Term 2 



1=1- 



o S T " 2J 



• s^ = i. 



For the second product, ( A13 ) gives 



T 



1 

T 



(A13) 

(A14) 

(A15) 
(A16) 

(A17) 



because only one of the T terms in the target has the angular momentum component 
in position i2j+i, as in the first factor. The result is 

1 



D (i) 



LoS 



T-1-2J 



LoS 



T-l 



T 



(A18) 
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Appendix A. 3. Canonical form of an even number of derivatives of Sq 

Here we discuss the derivation of an expression for the (2p)-derivative of Sq in which 
all terms are of rank T. We will carry out the derivation explicitly only for the terms 
proportional to Sq , since only these survive the contraction in (23). We can illustrate 
the problem appearing in the discussion of Section ^ following (25), with this example: 
take, say, T = 4 and p = 2 in (30a). In the sine-cosine representation, this is 

d 4 r 

(sin 4 fix) = fi 4 (32 sin 4 fix - 192 sin 2 Six cos 2 fix + 32 cos 4 fix) 



dx 4 



sin fix — 8 cos fix) . (A19) 



The first parenthesis contains products of "rank" 4 while the second contains products 
of "rank" 2. Direct differentiation, however, gives 



dx 4 



(sin 4 fix) = fi 4 [40 sin 4 fix- 192 sin 2 fix cos 2 fix+24 cos 4 fix] , (A20) 



where all the terms are of "rank" 4. The expansion of sin T fix in terms of sines 
and cosines of multip les a ngles before differentiation produced this effect. However, 
equations (A19) and ( |A2C| ) are identical: in (A19), multiplying the second parenthesis 
by cos 2 fix + sin 2 fix = 1 yields 8 sin 4 fix — 8 cos 4 fix. 

We need to perform the same regrouping of terms on (30 a, 306). This is very 
difficult to do in generic form. However, we should remember that these derivatives 
are to be multiplied by a factor proportional to Sq, which has components only in 
the radial direction. Consequently, all terms in the derivatives possessing tangential 
components (terms with cosines in the sine-cosine representation) will have no effect. 
The only surviving term in (30a, 306) is the one that has only radial components (all- 
sines term). Explicitly, this condition is T — 2k — 2q = in (30a) and T— 2k — 2q— 1 = 
in (306). 

One can demonstrate that in the regrouping process, the coefficient of the final 
all-sines term is obtained by simply adding the coefficients of all the all-sines terms of 
different rank T — 2k in the original sum (the same conclusion holds for the all-cosines 
term). The demonstration generalizes the example above, and involves multiplying 
the lowest ranked terms by cos 2 fix + sin 2 fix, combining with the next higher ranked 
terms, and repeating the process until all terms are the same rank. 

In (30a, 306), we thus need to add the coefficients of all S T 9 _ 2fe _ 2g , for T—2k—2q = 

and all k, in (30a), and the coefficients of all ^x-2k-2q-ii for T — 2fc — 2g — 1 = 
and all k, in (306), to get the coefficient of Sq. 

For (30a), the condition 2q — T — 2k, eliminates all the terms in the sum over q 
except the last one, so the coefficient of S P_2fc is 



(-l)L*J- fc Q (T-2fc)^(-l)L*J^ Q; JJf) = Q (T-2k)^. (A21) 

For (306), the condition 2q + 1 = T — 2k, also eliminates all the terms in the sum over 
q except the last one, so the coefficient of S^~ 2fe is 

( _ 1) m+^Q( T -2fc)^(-i)L^J^^-^ = Q {T _ 2k) 2 P (A22) 



18 



The result is the same for T even or odd. Thus, the term in the (2p) derivative of Sq 
that will survive the product by S is: 



(2p)cT 
a 



all-sines 



(-i)p n 2 p 

2 T ~ 1 



E 

fe=0 



(T - 2k) 



2p 



S^ = A(T,p) Si 



(A23) 



Appendix B. Mass functions fi(t]) 

In a general coordinate system, the expressions for the leading order of the radiative 
multipole moments are simply 



I«> = {rmxYxl 2 ■ ■ ■ x l { + m 2 44 2 • • • x\ l } gTF , 
J <;> = {e liab [m Y x a Y v\x^ ■ --x\ + m 2 x a 2 v h 2 x% ■ ■ -x\ l ] } 



STF 



(B24) 



Changing to the CM-coordinate system defined by the Newtonian relations (consistent 
with our bare- multipole approach) miXi + m2X 2 = 0, and with x = xi — x 2 , and 



m 2 
Xi = x, 

m 

TO2 

vi = v, 

m 



x 2 = 
v 2 = 



mi 

m 
mi 

m 



(B25) 
(B26) 



the moments become 
j_<i> 

i _ 



mi (^) +(-iym 2 m 

(m 2 \ l+1 , , /mi\ 
m M +(-l) J+i m 2 — 



x STF , 



i+1 



yU _1 [L o : 



.<«-!>! 



where L = /jx x v is the angular momentum vector, and x < ' > = x ll x 12 . . .x tl . We 
define p such that 



m 2 = p mi , m = (1 + p) mi , /U = 



(1 + P) 2 ' m ' 77 (1 + P) 2 ' 
with the convention m 2 < mi. Inverting the last expression, we have 
1 



2/7 



1 — 2?7— v/1 - 4?7 



(B28) 



(B29) 



In terms of and p, (B27) becomes 

f<'> - /; ^ 1 + x<'> = /; AM x <Z> 

r ~M (1 + pV-l STF — M M'/J X STF> 



[Lo 



r<i-l>l 



J STF 



[L°x 



<;-i>l 



I STF ' 



(B30) 



Notice that we have not used any information about the trajectory of the bodies in 
this section. This expression for //(??) is then valid for any system of two bodies. 

Figure 3 shows the dependence of fi on the reduced mass m for odd and even values 
of I. Limiting values for /; in the range r\ = (0, 0.25], are |/z(0)| = 1, (0.25) = 0, if I 
is odd, and /j(0.25) = l/2^ 2 , if I is even. 
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Figure captions 



Figure Al. STF multipolc moments displayed in an array with multipolc order 
increasing to the right, and the order of post-Newtonian corrections of each multipole 
moment increasing downward. Tail terms first appear at 0(e 3 ' 2 ) for each moment. 
For each I, ~£ <l> and f <l ~ 1> are grouped together. The polygons below indicate 
the moments and PN accuracies required to calculate the indicated quantity through 
2PN order. 



Figure A2. Fractional difference between series truncated at the labelled PN order 
and the exact series, as a function of v, for r) = (solid lines) and r\ = 0.25 (dashed 
lines). Innermost stable orbit for test body motion is at v Rj 0.4. In the equal-mass 
case, vanishing of odd-parity moments leads to degeneracy between adjacent PN 
approximations . 



Figure A3. The function |/;(??)|, for odd (dashed lines) and even (solid lines) values 
of I. Note that f 2 (r)) = 1. 



Figure A4. A comparison between the bare-multipole series (upper dashed curve), 
the "exact" numerical results from black- hole perturbation theory (solid curve), and 
the true PN series truncated to 3.5PN order (lower dashed curve). The three curves 
are plots of E/Eq as a function of orbital velocity v, in the physically relevant interval 
< v < 0.4. The comparison is valid in the test-body limit t)<l. 
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